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The quark number susceptibility, associated with the con- 
served quark number density, is closely related to the baryon 
and charge fluctuations in the quark-gluon plasma, which 
might serve as signature for the quark-gluon plasma formation 
in ultrarelativistic heavy-ion collisions. Beside QCD lattice 
simulations, the quark number susceptibility has been calcu- 
lated recently using a resummed perturbation theory (Hard 
Thermal Loop resummation) . In the present work we show, 
based on general arguments, that the computation of this 
quantity neglecting Hard Thermal Loop vertices contradicts 
the Ward identity and violates the thermodynamic sum rule 
following from the quark number conservation. We further 
show that the Hard Thermal Loop perturbation theory is con- 
sistent with the thermodynamic sum rule. 

Dynamical properties of a many particle system can 
be investigated by employing an external probe, which 
disturbs the system only slightly in its equilibrium state, 
and by measuring the response of the system to this ex- 
ternal perturbation. A large number of experiments be- 
long to this category [1] such as studies of various line- 
shapes, acoustic attenuation, and transport behavior. In 
all these experiments, one probes the dynamical behavior 
of the spontaneous fluctuations in the equilibrium state. 
In general, spontaneous fluctuations are related to cor- 
relation functions. They reflect the symmetries of the 
system, which provide important inputs for quantitative 
calculations of complicated many-body systems [2] . 

Recently, screening and fluctuations of conserved quan- 
tities, such as baryon number and electric charge, have 
been proposed as a probe of the quark-gluon plasma 
(QGP) formation in ultrarelativistic heavy-ion colli- 
sions [3-10]. These fluctuations are closely related to the 
quark number susceptibility (QNS) [3,7], which measures 
the response of the number density to an infinitesimal 
change of the quark chemical potential. 

The QNS has been considered within lattice gauge the- 
ory [11-13], in perturbative QCD [14], as well as effec- 
tive models of QCD [3,7-10]. Recently it has also been 
computed within the 2-loop approximately self-consistent 
$— derivable Hard Thermal Loop (HTL) resummation in 
Ref. [15] and within one-loop HTL perturbation theory 
(HTLpt) in Ref. [16]. Whereas the result of Ref. [16] is 
in good agreement with lattice results, Ref. [15] overesti- 
mates the lattice data. Note, however, that the entropy 
and pressure obtained within the 2-loop approximately 



self-consistent $— derivable method agree well with lat- 
tice data [17]. 

In Ref. [18] possible explanations for the different re- 
sults of the various HTL approaches have been discussed. 
As we will argue here, the neglect of the HTL vertices vi- 
olates the Ward identities and the thermodynamic sum 
rule following from the conservation of the quark num- 
ber [1,3,7]. We show, based on general arguments, that 
the HTLpt, on the other hand, is consistent with the 
thermodynamic sum rule. 

The HTL resummation technique allows a system- 
atic gauge invariant treatment of gauge theories at finite 
temperature and chemical potential taking into account 
medium effects such as Debye screening, effective quark 
masses, and Landau damping [19]. The generating func- 
tional which generates the HTL Green functions between 
a quark pair and any number of gauge bosons can be 
written [19,20] as 
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where K 11 is a light like four-vector, is the covariant 
derivative, m q is thermal quark mass, and ( ) is the aver- 
age over all possible directions over loop momenta. This 
functional is gauge symmetric and nonlocal and leads to 
the following Dirac equation 
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where we have suppressed the color index. In the HTL 
approximation the 2-point function, X ~ gT, (quark 
self-energy) is of the same order as the tree level one, 
Sq 1 ^) ~ Jf.~ gT (in the weak coupling limit g -C 1), if 
the external momenta are soft, i.e. of the order gT. The 
3-point function, i.e., the effective quark-gauge boson 
vertex, is given by gT M = 5(7^ + 5T^), where ST^ is the 
HTL correction. The 4-point function, g^^, does not 
exist at the tree level and only appears within the HTL 
approximation [19]. These JV-point functions, which are 
complicated functions of momenta and energies, are in- 
terrelated by Ward identities. This has important conse- 
quences for the QNS (and other quantities) if one aims 
at a consistent perturbative expansion. 

The HTLpt [21] is an extension of the screened pertur- 
bation theory [22] in such a way that it amounts to a re- 
organization of the usual perturbation theory by adding 
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and subtracting the nonlocal HTL term (1) in the ac- 
tion [19,23]. The added term together with the original 
QCD action is treated nonperturbatively as zeroth or- 
der whereas the subtracted one as perturbation. Equiva- 
lently, for calculating physical quantities this amounts to 
replace the bare propagators and vertices by rcsummcd 
HTL Green functions [19]. In this way, interesting quan- 
tities of the QGP have been computed [21,24-31] within 
the HTLpt. However, HTLpt suffers from the fact that 
one uses HTL resummed Green functions also for hard 
momenta of the order T. A systematic description of 
physical quantities requires an explicit separation of hard 
(~ T) and soft (~ gT) scales [19]. Following this ap- 
proach, a number of relevant physical quantities, such 
as signatures of the QGP, has consistently been calcu- 
lated [32]. A convincing reason for studying the QNS in 
HTL approximation is the fact that this approach has 
also been used for the free energy of the QGP, where 
different implementations of the HTL method and their 
validity have been discussed. The QNS may serve as an- 
other test for these approaches. 

The usual definition of the QNS for a given flavor which 
measures the response of the quark number density p q to 
an infinitesimal change in the quark chemical potential 
H q + S{i q is given by [1,3,7]: 
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Since the quark (baryon) number is a conserved quantity, 
the QNS obeys a Kubo relation [1] 
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where jo is the time component of the quark current = 
i/>7 M V- The equivalence of the static susceptibility (3) and 
the static limit of the dynamical susceptibility (4) follows 
from [1] 
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with the quark number operator Af q — J d 3 xjo(0, x) and 
volume V . The derivative with respect to /x g in (5) ac- 
cording to (3) agrees with (4) only if the number opera- 
tor is a conserved quantity and hence commutes with the 
HamiltonianH [1]. 

Introducing the vector meson correlator S^i,(t,x) = 
{jii{t,x)jv(0,0)) and taking the Fourier transform of 
500(0, x), the QNS according to (4) can be written as [33] 
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This equation is known as a thermodynamic sum rule as it 
relates the thermodynamic derivatives to the correlation 



function due to the symmetry of the system. Examples 
of such sum rules for various systems can be found in 
Ref. [1]. The fluctuation - dissipation theorem relates 
Soo to the vector meson self-energy n o via [1] 
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An example for this self-energy within the 1-loop order 
with effective Green functions is shown in Fig.l, where 
the blobs indicate effective quark propagators and quark- 
meson vertices. 

The time independence of Af q , corresponding to the 
continuity equation 9 M j p = 0, implies that the quantity 
J d 3 xS o(0,x) in (4) is time independent. This is guar- 
anteed if Soo (u), 0) cx S(ui) [1]. 
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Fig.l: The 1-loop vector meson self-energy diagrams. 

In usual perturbation theory the equivalence of (3) and 
(4) follows from the fact that the loop expansion agrees 
with an expansion in the coupling constant. For example, 
to lowest order the number density follows from the tad- 
pole diagram (containing only a bare propagator and ver- 
tex) at finite temperature and chemical potential, which 
corresponds to the 1-loop polarization tensor of the vec- 
tor correlator (7). For the HTL rcsummation the loop 
expansion and the coupling expansion are not equivalent 
as the first one includes higher order effects through re- 
summation. This mixing can cause an incompleteness in 
physical quantities computed in a given loop order, as we 
shall discuss below. 

Now we would like to discuss how the QNS can be cal- 
culated consistently within the HTL approximation. Our 
starting point is the most general expression of the quark 
self-energy [34], which enters the full quark propagator. 
In the rest frame of a medium in the chiral limit it reads 
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where the scalar quantities a and b are functions of the 
energy ko and of the magnitude of the spatial momentum, 
k, and K = (k ,k). The full propagator is given by 



Sp 1 (K) = If - S (K) = (1 + a) &7 
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Using the hclicity representation [29] it can also be writ- 
ten as 
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Charge invariance demands D±(—k ,k) — —D^(ko,k), 
which implies a(— ko,k) = a(ko,k) and b(— ko,k) = 
— b(ko, k). The zeros of D± describe the in-medium prop- 
agation or quasiparticle (QP) dispersion relation of a par- 
ticle excitation with energy uJ+{k) and of a mode called 
plasmino with energy, uj-(k), which exists only in the 
medium [29,35]. In addition, D± contains a discontinu- 
ous part corresponding to Landau damping (LD) if the 
quark self-energy has an imaginary part. 

Using these general properties of the quark propagator, 
the spectral representation of the in-medium propagator 
(10) can be written as 

p±(u, k) = — — Im— J— = 1Z±(ui, k)S(uj — u>±) 

+K T (-uJ 7 k)5(u J + u JT ) + pf s (u} 7 k) , (11) 

where 1Z± = (dD±/doj)~ 1 are the residues correspond- 
ing to the poles of the propagators, and p± IS are discon- 
tinuities. In the HTL-approximation the discontinuities, 
given explicitly e.g. in Ref. [16], are located below the 
light cone (k 2 < k 2 ), where Landau damping occurs. We 
also note that for each pole the spectral function will in- 
volve a ^-function. The Ward identity relates the quark 
propagator to the quark-meson vertex via 



(K 2 - K^T^-K^K^Ki - Ki) 
= S F \K 2 )-S F 1 (K 1 ). 
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In the HTL approximation the effective vertices are dis- 
cussed in Ref. [16]. 

Using the effective propagators and vertices given in 
(10) and (12), the contribution of the first diagram in 
Fig. 1 can be written as 

Imn; (w, 0) = 4N f N c n (l - eP") J J dx J dy 

x 5(u-x-y)n F (x) n F (y) |(1 - f(x, y)f 
x P+{x, k) p_ (y,k) + -\- [Im F+ (x, k) 

U> 1 7T 



where F±(z, k) = D±(z, k) + z^f k, z is the energy of the 
internal quarks, f(x,y) = {x + y)/uj, Nf is the number of 
quark flavors, and N c is the number of colors. The contri- 
bution of the tadpole diagram, ImnQ (ui, 0), in Fig. 1 is 
identical apart from an opposite sign of the second term 
inside the curly braces in (13). Using the explicit form 
of the spectral functions (11), the total contribution of 
Fig. 1 is 

rmUoo (w,0) = [Imng (w) + Imll£ (w)] = 4N f N c ir 
(1-^)%) / ^[Ul(u + ,k)n F (u + )n F (-u; + ) 

J (27T) 

+ H 2 _(u-,k)n F (u-)n F (-u-)] , (14) 



where the LD contributions explicitly vanish. Imiloo in 
(14) contains only the QP contributions following from 
the dispersion relations, D± = 0, and the requirement 
for the conservation law, Soo oc S(uj) 7 is satisfied. This 
necessitates the use of effective 3- and 4-point functions 
and also the modification of the charge operator, as 
"0(7o + STq)-)P, by inclusion of the effective 3-point func- 
tion, in contrast to Ref. [15,18]. This non-local modifi- 
cation makes the charge operator different from that of 
(3) and (4) but is essential to satisfy the Ward identity 
which is related to the gauge invariance and conservation 
law. This implies that the number density has to be cal- 
culated with propagators and vertices related by Ward 
identities at finite T and \x q . 
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Fig. 2: The QNS in different Ij^Teximations (see text) 
divided by the lowest order result (free QNS) as a func- 
tion of T/T c with TJAj^ = 0.49 [36]. In each band the 
lower curve corresponds to the choice of the renormaliza- 
tion scale p, = 2nT and the upper one to p, = A-kT. The 
squares represent recent lattice data [13]. 



x p_(y,k)+ImF_(y,k)p+(x,k) ]}, (13) (6) 



Now the 1-loop QNS follows from combining (14) and 
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which is identical with the one-loop HTLpt result [16], if 
the residues are obtained from the HTL resummed prop- 
agators. In this case the functions 1Z± in (15) are given 

by 



i4(fc)-fc 2 

*±(»±> k ) = 2m 2 
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with the HTL quark dispersion relations u)±(k) and the 
HTL effective quark mass m q = g(T)T/\^6. The re- 
sult (band between the solid line denoted by QP in 



3 



Fig. 2) agrees well with recent lattice data [13]. (Here 
we used the 2-loop result for the running coupling con- 
stant g(T) in contrast to Ref. [16], where the 1-loop ex- 
pression was adopted. For the renormalization scale ft, 
entering the running coupling constant, we choose two 
different values leading to the bands in Fig. 2.) How- 
ever, it overcounts the leading order perturbative re- 
sult (dotted band denoted by PLO in Fig. 2) given as 
Xq/Xf = (1 — 2a s /-K + • • •) [18]. We note that this can 
be cured by going to the calculation at 2-loop order. Re- 
cently, it has been shown that the 2-loop order calcula- 
tion [27,28] for the thermodynamic potential reveals the 
correct inclusion of the leading order effects but happens 
to be very close to the 1-loop result. 

The physical relevance of the plasmino branch has been 
discussed in connection with the dilepton production in 
rclativistic heavy-ion collisions [29,35]. Using the HTL 
approach its contribution to the QNS, i.e. the second 
term of (15), amounts only to less than one percent of 
the total susceptibility. This is caused by the fact that 
within the HTL approximation the spectral strength of 
the plasmino branch vanishes exponentially for hard mo- 
menta, k » gT, and that we integrate in (15) over all 
momenta. 

Now, employing only resummed propagators and bare 
vertices as in Rcf. [15], the QNS is obtained from (13) 
by setting the vertex correction functions, f(x, y) and 
F±(x, y), equal to zero. Then the QNS contains LD con- 
tributions coming from the discontinuity of the 2-point 
HTL function in addition to the QP contribution in (15). 
The result is shown in Fig. 2 by the dashed band de- 
noted by QP+LD which clearly overestimates the lat- 
tice data. Although the charge operator is local in this 
approach [18], 5oo(^,0) is not proportional to S(w), im- 
plying that dp q /dt 7^ 0. This amounts to a violation of 
the Ward identity, i.e, gauge invariance, for the effective 
HTL Green functions and also sacrifices the conservation 
law and hence the thermodynamic sum rule. Note, how- 
ever, that it has been argued that the gauge dependence 
of the 2-loop effective action within the ^-derivable ap- 
proximation affects only higher order contributions [37] . 

In summary, we have discussed that the QNS can be 
computed consistently within the HTLpt in order to sat- 
isfy the requirement of the thermodynamic sum rule. 
This requires to consider 3-point and 4-point HTL func- 
tions, related to the HTL propagator by Ward identities. 
As discussed above, the QNS within the one-loop order 
overestimates the PLO, which, however, can be cured 
by extending the calculation to 2-loop orders in HTLpt. 
This would, indeed, be a challenging as well as nontrivial 
task due to its complexity. However, if it turns out to be 
a small correction, e. g., like in the case of the quark-gluon 
thermodynamic potential, the 1-loop HTLpt result pre- 
sented here, agreeing well with lattice data, would be a 
useful approximation. 
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